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ABSTRACT 

We study the feedback between heating and cooling of the intra-cluster 
medium (ICM) in cooling flow (CF) galaxies and clusters. We adopt the popular 
view that the heating is due to an active galactic nucleus (AGN), i.e. a central 
black hole accreting mass and launching jets and/or winds. We propose that 
the feedback occurs with the entire cool inner region (r < 5 — 30 kpc), where 
the non-linear over-dense blobs of gas with a density contrast p/p a > 2 cool fast 
and are removed from the ICM before experiencing the next major AGN heating 
event. We term this scenario cold-feedback. Some of these blobs cool and sink 
toward the central black hole, while others might form stars and cold molecular 
clouds. 

We derive the conditions under which the dense blobs formed by perturbations 
might cool to low temperatures (T < 10 4 K), and feed the black hole. The 
main conditions are found to be: (1) An over-dense blob must be prevented from 
reaching an equilibrium position in the ICM: therefore it has to cool fast, and the 
density profile of the ambient gas should be shallow; (2) Non-linear perturbations 
are required: they might have chiefly formed by previous AGN activity; (3) The 
cooling time of these non-linear perturbations should be short relative to few 
times the typical interval between successive AGN outbursts. (4) The blobs 
should be magnetically disconnected from their surroundings, in order not to be 
evaporated by thermal conduction. 
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Subject headings: galaxies: active - galaxies: clusters: cooling flows - galaxies: 
clusters: individual: M8 7/ Virgo 

1. Introduction 

The observations of Galaxy Clusters with the last generation X-ray satellites Chandra 
and XMM-Newton have shown the remarkable lack of the large amount of cool gas predicted 
by the old version (ten years ago) of the cooling flow (CF) model (e.g. Peterson et al. 2003; 
Tamura et al. 2001b; Molendi & Pizzolato 2001). The most straightforward explanation is 
that the intra-cluster medium (ICM) in CF clusters must be heated by some mechanism 
(e.g. Fabian 2003). Probably the most popular heating engine is an active galactic nucleus 
(AGN) residing at the center of the cluster dominant galaxy, often a cD (see e.g. Binney & 
Tabor 1995; Tucker & David 1997; Ciotti & Ostriker 2001; Binney 2004; for more references 
see Peterson, Kahn, Paerels, Kaastra, Tamura, Bleeker, Ferrigno, & Jernigan 2004 and the 
review by Begelman 2004; see Fujita &: Reiprich 2004 and Fabian, Voigt, & Morris 2002 for 
some problems connected with this scenario). Other heating mechanisms, like supernovae 
and heat conduction from the cluster outer regions were shown to be problematic for several 
reasons. It is sufficient to mention that we are looking for a unified mechanism to heat the 
ICM in CFs, from galaxies to large clusters. Heat conduction cannot work in galactic CFs 
because there are no large heat reservoirs, and probably there are not enough supernovae in 
elliptical galaxies to supply a sufficient amount of heat. 

The new X-ray observations show that the mass cooling rate to low temperatures is far 
below predictions by older versions of the cluster CF model (Fabian 1994), but is compatible 
with low- mass cooling rate models (see Binney 2004, on the expectation for that result). 
Among them, the moderate CF model (Soker et al. 2001; Soker Sz David 2003; Soker 2004) 
is different from many earlier proposed processes whose aim is to prevent CF in clusters of 
galaxies altogether. In the moderate CF model small but non-negligible quantities of gas are 
cooling to low (< 10 4 K) temperatures. 

Most models of AGN heating agree in that there is some sort of feedback between the 
heating and the radiative cooling. Many of these models result in intermittent AGN activity. 
There are two approaches to the feedback between the cooling ICM and the AGN. In the first, 
the ICM does not cool below X-ray emitting temperatures; the AGN accretion is Bondi-like, 
and is determined by the ICM properties very close to the central black hole (e.g. Churazov 
et al. 2002; Nulsen 2004; Omma & Binney 2004). Nulsen (2004) for example, studies the 
accretion of gas at the virial temperature. We term this type of models hot feedback. 
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In the second approach, the black hole accretes cold gas, but the mass cooling rate is 
much below that in old versions of the cooling flow model. In this case the feedback takes 
place within a region extending to a distance of ps 5 — 30 kpc from the cluster center. Such is 
the moderate cooling flow model (Soker et al. 2001; Soker 2004). We term this cold feedback 
model. The two types of feedback models imply other differences. 

1. In the hot feedback models the optical filaments observed in many CF-clusters (Heck- 
man et al. 1989), and the cooler molecular gas detected via CO observations (Edge 
2001; also Edge & Frayer 2003; Salome & Combes 2003; Salome & Combes 2004) come 
from stripping gas from cluster galaxies. In the cold feedback models the cold gas may 
come from the cooling ICM as well. 

2. In the cold feedback models some X-ray emission from gas at temperatures < 10 7 K is 
predicted, but at a level more than an order of magnitude below that in old versions 
of the cooling flow model (Pizzolato, Behar, & Soker in preparation). 

3. The feeding of the central black hole with cold gas in the cold feedback models makes 
the process similar in some aspects to that of AGN in spiral galaxies. 

4. It is possible, although not required, that in the cold feedback model the feedback does 
not only keep energy balance, but mass balance as well. Namely, a substantial fraction 
of gas that cooled to low temperature is injected back to the ICM and heats up as it 
is shocked (Soker & Pizzolato 2005). It is possible that during a fraction of the time, 
most of the cooling gas forms molecular clouds and stars, rather than being injected 
back to the ICM at high speeds. 

To set the stage for the discussion to follows, it will be constructive to present the 
temperature profiles of some CF clusters (Figure 1). In the presently proposed model a 
flat temperature profile in the inner region is expected, at least during some fraction of the 
duty cycle. It is that flat and low temperature profile which facilitates the formation of cold 
blobs. This is contrary to some hot feedback models, where the accreted gas comes from 
the immediate AGN neighborhoods. It is our opinion that some of these models may have 
problems in accounting for the temperature profiles presented in Figure 1. For example, 
the repeated heating by Omma &: Binney (2004) does not heat the outer region as much 
as the inner one. It seems as if the outer region will cool to low temperature, unlike the 
higher temperature of the outer regions in clusters. Ruszkowski Sz Begelman (2002) take 
both AGN heating and heat conduction. We find their temperature profile in the inner 
several kpc region to be too steep (see their Figure 1) compared with real clusters. In the 
temperature profile obtained by Hoeft & Briiggen (2004, their Figure 5) we see two problems. 
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First, contrary to their claim, we do not think that in the outer CF region their profile fits 
that of A2052 (taken from Blanton et al. 2003). Second, after a time of ~ 1.2 x 10 10 years, 
the temperature in the inner « 2 kpc drops to <C 10 7 K, contrary to the observations. There 
is a fine tuning problem, in that presently all clusters do not show this drop, but they will 
in a few 10 9 years. 

The considerations above, among other arguments, motivate us to consider the cold 
feedback model. In the present paper we examine the two types of feedback heating. By 
analyzing published models and by comparing them with the observations, we constrain the 
parameters space of the different models, and make some predictions which can be tested 
with future observations. Our proposed scenario is presented in § 2. In § 3 we calculate the 
evolution of a dense blob, assuming the dense blobs are magnetically disconnected from their 
environment, e.g., as in the magnetic flux loop model of Soker (2004) where heat conduction 
occurs only within the cold blob. Readers interested only in the basic scenario, results, and 
predictions, can skip § 3 and go directly from § 2 to our summary in § 4. 



2. The Proposed Cold Feedback Scenario 

The cold feedback scenario entails a cycle in the cooling/accretion activity. We suggest 
that this cycle starts with a major AGN outburst, which injects a huge amount of energy 
into the ICM. This event triggers the formation of a wealth of dense blobs. It is important to 
realize that these blobs are non-linear perturbations of the ICM, and may be distributed with 
a wide spectrum of densities. These blobs are denser than the surrounding medium, and fall 
to the black hole. If these blobs have an initial angular momentum, they do not freely fall, 
but may form an accretion disc. As we shall demonstrate (§ 3.1), for the expected angular 
momentum distribution this is not a relevant complication. The dense and cool blobs are 
the fastest to fall, and therefore are removed first from the ICM. The accretion history may 
be difficult to predict, since it depends on several factors. For most of the time the small 
blobs may accrete approximately in a steady-state, but some processes may intervene to 
modify this. For instance if most of the blobs are formed at the same distance to the AGN 
with a sharply peaked density spectrum, most of the blobs might accrete simultaneously. 
Moreover, an accretion disc may form and undergo some kind of instability. In either case, 
a sudden "catastrophic" accretion episode on the black hole is expected, resulting in a new 
AGN outburst, which restarts the cycle with a fresh injection of blobs. We examine here 
below some details of the suggested process. 

In this duty cycle some of the gas cools to low temperatures (< 10 4 K) before the next 
major heating, while the rest is heated back to a relatively high temperature. We differ from 
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many previous models in that in the moderate CF model a substantial fraction of the ICM 
gas cools to low temperatures. The presence of a detectable amount of gas cooling below X- 
ray emitting temperatures is a prediction of this model. Indeed, in the CF cluster Abell 2597 
both extreme-UV and X-ray observations indicate a mass cooling rate of ~ 100 M Q yr" 1 , 
which is ~ 0.2 of the value quoted in the past based on ROSAT X-ray observations (see the 
discussion in Morris & Fabian 2005). In the CF cluster Abell 2029, Clarke et al. (2004) find 
a substantial amount of gas at a temperature of ~ 10 6 K; a CF model gives a mass cooling 
rate of ~ 50 M & yr _1 . 

The AGN outburst interacts in a very complicated fashion with the ICM (e.g. Begelman 
2004), e.g. it heats and inflates radio bubbles, which rise buoyantly in the ICM. The ICM 
itself is displaced and thickened by the rising bubbles, as shown by the enhanced X-ray 
brightness (e.g. Blanton et al. 2001). A non- homogeneous thickening may result in the 
formation of a multi-phase gas inside or near the radio lobes, which harbor relatively strong 
magnetic fields, up to few tens \iG. 

The magnetic field inside the radio lobes is related to a fundamental issue for the cold 
feedback model, namely the efficiency of thermal conduction. A highly efficient conduction 
would evaporate the cold gas blobs before they can accrete on the AGN, and is therefore 
incompatible with the cold feedback model. 

It is well-known that magnetic fields are able to suppress thermal conduction, and that 
the degree of suppression strongly depends on their topology (see e.g. the discussion in 
Narayan & Medvedev 2001). 

For our purposes, we assume that there is essentially no heat conduction between the 
over-dense blobs and their surroundings, i.e., the over-dense blobs are magnetically discon- 
nected from their environment, as in the model of Soker (2004). Effectively, for our blobs of 
radius a ~ 10 — 100 pc, we find from Figure 3 of Nipoti & Binney (2004) that the effective 
heat conduction should be < 0.001 times the Spitzer (1956) value in order for the blobs not 
to be evaporated. 

This suppression factor is somewhat high, but consistent with some recent observations. 
Molendi (2002), e.g., finds evidence of a gas component cooler than the ambient gas inside 
the radio lobes of M87. The very existence of these cold pockets led Molendi (2002) to 
estimate a conduction suppression factor < 0.01 with respect to the nominal Spitzer (1956) 
value. 

The conclusion is that the radio lobes, permeated by magnetic fields, may be "safe 
corridors" where the blobs may accrete on the AGN without being evaporated by conduction. 
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The coexistence of gas phases at different temperatures is not exclusive of M87. In their 
analysis of the group NGC 5044, Buote et al. (2003) find evidence for a moderate multi- 
phaseness. Their data are well fit by a two temperature model, a relatively cool component 
with T coo i pa 0.7 keV and a hot component at T\ lo t ~ 1.4 keV. These temperatures seem to 
coexist in the inner pa 30 kpc, with the cooler component dominating in the inner pa 10 kpc. 
The cool component has a sizeable filling factor: /# pa 0.5 for r < 25 kpc and fs ~ 0.1 at 
r > 30 kpc. 

We therefore consider the inner region, r < 5 — 30 kpc of cooling-flow clusters to posses 
an ICM with non-linear perturbations, i.e. dense blobs spread within it. The fate of a blob 
depends on the relative magnitudes of its cooling time t coo i> the time interval to the next 
AGN heating event, and on the time tf a n the blob takes to accrete on the central black 
hole. The time scale tf a n depends on several factors. The first is the blob's over-density with 
respect to the ambient gas. The second is the blob's angular momentum: a high angular 
momentum prevents the blob from approaching the black hole altogether. We shall discuss 
the issue of angular momentum more fully in § 3.1: for the time being we assume that the 
content of angular momentum of a blob is small. Under this hypothesis, the relation between 
the time scales t coo i and tf a n mainly depends on the blob's density p relative to that of the 
ambient medium p a . A sinking blob has t coo \ < tf a n, with the extreme case of t coo i <C £f a u, 
where the blob significantly cools down on a very short time scale, and almost free-falls to 
the center. If tf a u < t coo i the blob would not cool much, and will reach a position where it 
is as dense as the ambient gas; if t coo \ ~ £f a n the blob sinks fairly slowly, owing to its small 
over-density. 

In this way, the ICM efficiently disposes of the cooler phases, i.e., the highly non-linear 
over-dense perturbations. 

If a blob cools isobarically in pressure equilibrium with its surrounding, its cooling time 
scales as /-cool 

oc p . Later on, when the blob temperature drops below T pa 0.1 keV, 
the sound waves become too slow to keep up with the outer pressure, and cooling occurs 
isochorically: t coo i oc p -1 (e.g. Burkert & Lin 2000). In either case the denser — and the 
cooler — is a perturbation, the faster it cools, and the more efficiently it is removed from 
the ICM. 

The discussion above leads us to the proposed scenario. Non-linear over-dense blobs of 
gas, Sp/p a > 2, i.e. p/p a > 3, cool on short time scales such that they are removed from the 
ICM before the next major AGN heating event. Some of these blobs cool and sink toward 
the central black hole. Other non-linear perturbations may form stars, as is inferred in some 
CF clusters, e.g., in the CF cluster A1068 the cooling rate within r pa 30 kpc is about equal 
to the star formation rate there (Wise et al. 2004; McNamara et al. 2004). The dense blobs 
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that sink to the center feed the AGN. The feedback is with the entire cool inner region, and 
not only with the gas close to the black hole. Any over-cooling taking place in the inner 
region, where the temperature profile is flat, will lead to many small and dense blobs, which 
feed the AGN. 

We assume that the cold blobs are magnetically disconnected from the surrounding, 
so that most of them can survive long enough to be delivered to the central AGN. Indeed, 
while some blobs might certainly be evaporated, yet the inability of thermal conduction to 
re-heat all the cold gas is testified by the presence of sizeable amounts of molecular gas in the 
central few kpc of some clusters (NGC 1275/Perseus: Inoue et al. 1996). In a recent paper 
Wilman et al. (2005) resolved a ring of molecular gas with radius of 50 pc from the center 
of NGC 1275 at the center of the Perseus cluster. Adding the presence of molecular gas at 
distances of up to ~ 10 kpc, their finding shows that cold gas originating at large distance 
in the cluster, but still within the low temperature region, can feed the central black hole. 

This ends our simple demonstration that highly over-dense blobs can be accreted to the 
central black hole before the next major AGN heating event. More detailed calculations are 
presented in the next Section; readers not interested in them may skip directly to the last 
Section. 



3. Nonlinear Evolution of the Blobs 

In this Section we consider the evolution of a single blob. Let V and S be its volume 
and cross-section, respectively, and p its mass density. The blob is subjected to the overall 
gravitational acceleration g, to the hydrostatic buoyancy force and to the drag force. Its 
equation of motion may be written (e.g. Loewenstein 1989; Kaiser 2003) 

pV — = g {p- p a )V - — S p a vv, (1) 

where p a is the mass density of the ambient gas and C is the dimensionless drag coefficient. 
For C we shall adopt the value C ~ 0.75 used by Kaiser (2003) and derived from the 
numerical simulations by Churazov et al. (2001). It will be convenient to express Equation (1) 
in terms of the blob over-density 5 with respect to the surrounding medium 

5= P -^. (2) 

Pa 

Eliminating p from Equation (1), we may rewrite it as 

dv 5 3C v v . . 
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where we have calculated V and S assuming that the blob is a sphere of radius a. If a, 5 
and g are constants, the last equation is easily integrated: 

JF + i 

where 



e ■ 
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is the characteristic time in which the blob attains its terminal velocity 
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Note that if the gravitational acceleration g is relatively small, then a blob will take more 
time t oc (7 -1 ^ 2 to accelerate to a smaller terminal velocity Vt oc (? 1 ^ 2 . This explains why in 
Figure 2 the terminal-velocity fall time drops below the free-fall time at small radii, where 
g is small; the blobs starting there take a comparatively long time to accelerate. 

In order to calculate the evolution of a blob, we need three further equations. The first 
equation simply relates the position and the velocity of the blob: 

dr , . 

The second equation is the blob mass conservation 

p a 3 = constant. (8) 

The last equation provides the blob energy balance. If E is the blob total internal energy 
and P is its pressure, the combined first and second law of thermodynamics 

P dS = -n e n H AVdt (9) 

yield the variation of the blob entropy 

s = h* (-41 (io) 



2 V™ 5/3 , 

under the radiative loss of the amount of energy dQ = n e nu A V dt in the time interval dt. 
This loss depends on the cooling function A(T) and on the electron and hydrogen densities 
n e and nu of the gas inside the blob. We now suppose that the blob is instant by instant in 
pressure equilibrium with the ambient gas, whose pressure is P a : 



P = P a - 



(11) 
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This assumption holds until the blob sound crossing time is short in comparison to the blob 
cooling time, and therefore the acoustic waves can efficiently pressurize the blob (see e.g. 
Burkert & Lin 2000). We shall have to check a posteriori the validity of this hypothesis. We 
further assume that the ambient gas is in hydrostatic equilibrium. With these hypotheses, 
the pressure variation dP experienced by the blob is only due to its drift across the cluster 
gravitational potential $: 

dP = dP a = -p a d$ = -p a dtv-V<S> (12) 

With the help of this equation we differentiate Equation (10) and plug the result into Equa- 
tion (9): after some straightforward algebra we obtain 

dlnn n a g ■ v 2 n e n# A(T) 
dt n c 2 s 5 n kT ' 

where c 2 s = 5 kT/3 pm p is the (squared) adiabatic sound speed inside the blob, k is the 
Boltzmann constant, g = — V$ is the gravitational acceleration and n a = p a /pm p is the 
ambient total number particle density. 

Using the definition 2, we may rewrite Equation (13) in terms of the over-density 5: 
d8 g-v 2 ^ „, A(T) M 

-j t = + ^PePH n a (1 + 5) 2 - (1 + 5) v ■ V lnn a , (14) 

where the coefficients p e and pu are defined by p e = n e /n and pu = uh/ti. The first term 
on the right-hand side of Equation (14) is the adiabatic compression owing to gravity; the 
second term describes the blob thickening at constant pressure as its temperature cools on 
account of the radiative losses; the last term is proportional to the slope on the background 
ambient gas, and is the reduction of the density contrast owing to the denser layers the 
blob sinks through. This last term is important in what it may hamper the growth of large 
over-densities in a strongly stratified ambient gas. 

Equations (2, 3, 7, 8) and (14) provide a closed system of equations for the blob evolution, 
once we have specified the distribution of the ambient density p a , the cluster gravitational 
potential $ and a suitable set of initial conditions. 

It is important to remark that in our derivation we have assumed that the ambient 
gas behaves as a static background without cooling. This approximation is justified as far 
as the blobs cool faster than the ambient gas. This assumption certainly holds if the blob 
remains over-dense, but may break down if \8\ 1 (Loewenstein 1989). Therefore, our 
results concerning blobs only slightly over-dense with respect to their surroundings require 
some care. 



(13) 
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We apply the above equations to calculate the evolution of a blob in the cool core cluster 
M87. We choose this cluster because we know to a good approximation its gravitational 
mass, its temperature and density structure; besides, the central AGN here shows a complex 
interaction with the surrounding medium (e.g. Belsole et al. 2001; Molendi 2002). Where 
necessary, however, we shall extend our considerations beyond the particular features of 
M87. For the ambient density p a , temperature T a and gravitational acceleration g of M87 
we assume the functional forms provided by Ghizzardi et al. (2004). We have employed the 
cooling function A(T) given by Sutherland & Dopita (1993) for a solar-abundance thermal 
plasma. As far as the initial conditions, we assume that the blob has an initial over-density 
5q > 0, and is released from rest at the distance Tq from the cluster center. We have solved 
numerically the system of Equations( 2, 3, 8, 14, 7) and these initial conditions with a 
step-adaptive fifth-order Runge-Kutta algorithm (Press et al. 1992). 

Figure 3 plots the over-density evolution of some blobs, which only differ in their initial 
over-densities So', all of them have initial size ao = 100 pc, and have been released from 
rest at tq = 20 kpc from the cluster center. It is apparent that the qualitative behavior of 
the solution critically depends on 5q: if 5q is below a critical threshold 5c the over-density 
decreases, else if 5 > 5c it diverges. A blob with 5 < 5c soon evolves to \5\ <C 1, becoming 
hardly distinguishable from the ambient background gas, whose cooling and bulk velocity 
have been neglected in our approximation. Therefore, our result concerning these blobs may 
be looked at suspiciously, and for consistency we conservatively disregard their contribution 
to the feeding of the central black hole. These blobs are most likely to get pinned by the drag 
force to the average bulk motion of the ICM, and their successive evolution would require a 
more detailed analysis (see e.g. Nulsen 1986; Loewenstein 1989). To our concern, since the 
background medium has been assumed static, we may assume that such blobs have been 
thermally stabilized. 

The existence of a critical over-density 5c is not difficult to understand (Balbus & Soker 
1989; Loewenstein 1989; Tribble 1991). As a blob sinks, it crosses thicker and thicker layers 
of ambient gas, and its density contrast 5 with respect to its surroundings would reduce. If 
this effect is not contrasted, the blob would stop its fall at an equilibrium distance from the 
center where its density is the same as the ambient gas, i.e. 5 = 0. The only way to overcome 
this effect is to cool fast. Since this process is isobaric, the blob density is enhanced, and the 
blob might reach the cluster center before its over-density has significantly reduced. Since 
the blob cools on a time-scale i coo i, a blob condenses and reaches the center if t coo i < tf a u. 
The equality between these time-scales defines the critical value 5c for the blob over-density. 
This explains the qualitative behavior shown in Figure 3. 

Figure 4 plots the evolution of the position and velocity of two blobs, one slightly over- 



- 11 - 



critical (5 > 5c) and one slightly under-critical (5 < 5c)- In both cases, the blob is initially 
in free-fall, but after a transition time of the order of few r -given by Equation (5)- the drag 
limits the velocity to the terminal value v t given by Equation (6). An under-critical blob 
attains with non-zero velocity its equilibrium position (where 5 = 0), which is overshot. The 
blob starts a series to oscillations about this position, which are quickly damped by the drag 
force (e.g. Balbus & Soker 1989; Loewenstein 1989; Tribble 1991). Eventually, the blob cools 
and becomes pinned to the bulk of the ambient gas. An over-critical blob is always denser 
than the surrounding gas, the drag plays a minor role, and the blob nearly free falls all the 
way down to the center. 

The evolution of the blob temperature is presented in Figure 5. The under-critical 
blob has an initial radiative cooling time of t coo i ~ 1.8 x 10 8 yr, and before this time the 
blob temperature does not change much. For such a blob the gravitational heating is more 
important than the radiative losses: the blob is heated up to the local ambient temperature, 
and becomes thermally stable. In a denser blob, on the contrary, the gravitational heating 
is insufficient to overcome the radiative cooling. After a time t coo i ~ 1-2 x 10 8 yr, therefore, 
the blob cools very fast and almost free falls to the cluster center. 

The initial release radius is important for the nonlinear development of a blob. Figure 6 
plots the over-density evolution for some blobs with the same intrinsic properties, but released 
at different radii. It is apparent that the farthest blobs never develop large over-densities. 
This fact is related to the density profile of the ambient gas (shown in Figure 7). For a 
fixed over-density, a far blob has a relatively small density -in absolute terms- and as it 
sinks it soon gets embedded in the denser ambient gas: the over-density reduces and the 
blob is thermally stabilized. The situation is different if the blob has been released at a 
small distance to the center, where the ambient density profile is flat (Figure 7). As in 
our discussion following Equation (14), in this case as the blob's over-density growth is not 
hindered by the ambient gas gradient, and even a moderate initial over-density may evolve 
towards larger values. It is worth to stress the importance of this density stratification effect, 
as it relates the size of the central density plateau to the amount of allowed cold feedback in 
our model. Indeed, in the density plateau essentially all the over- densities may evolve to the 
non-linear regime, cool down and accrete on the AGN to provide the heating feedback. On 
the other hand, only the largest over-densities born outside the central plateau may evolve 
this way. In the case of M87 this plateau extends out only to r f=s 5 kpc (Figure 7): more 
typical cooling flow clusters, like A2052, have r rs 30 kpc (Blanton et al. 2001); this allows 
the feedback between the AGN and a larger fraction of the cluster's gas. 

Following Tribble (1991) we estimate the critical over-density 5c as a function of the 
other parameters. Equation (14) evaluated at t = with our initial conditions yields 5 > 0. 
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If the blob is stable, this initial trend must be reversed, and there must be an instant t = t 
for which 8(t = t) = 0. If the blob is unstable, on the other hand, 5 > all the way. 
Equation (14) evaluated at t — t may be written 



gr d In n a 



c? din r 



a 



t a v A(T a ) 

~ Mr)' [ j 



where c\ = 5 kT a / 3fim p and t a = 5 kT a /2 fi e fi H n a A(T a ) are respectively the sound speed and 
the isobaric cooling time of the ambient gas; in pressure equilibrium the blob temperature 
T and the ambient gas temperature T a are related by T = T a /(1 + 5). We must evaluate 
the right-hand side of Equation (15) at t — t . If r <C t (where r is given by Equation 5), we 
may substitute v with its terminal value given by Equation (6). If we suppose that the other 
quantities on the right-hand side of Equation (15) are not too different from their initial 
values, the last equation may be rewritten 



5 1 / 2 a \ 3C 



1/2 



gr dlnn 



a\V2 A(T a ) 



rJ A(T) ' 



(16) 



The value of 5 provided by this equation is the critical threshold 5c between a thermally 
stable and a thermally unstable blob. In the non-linear regime 5 ^> 1, by approximating 
the cooling function by a power law A oc T a , we obtain 5c oc a 1 ^ 3-2 ^, which coincides with 
the expression given by Tribble (1991) if a = 1/2, i.e. if thermal Bremsstrahlung is the 
chief coolant. Figure 8 plots 5c as a function of a for different release radii r (see also the 
Figure 1 of Tribble 1991). In the upper branch of the 5c — a plot we find 5c oc a 1 / 3 owing 
to the weak dependence of the cooling function on T in the temperature range typical of 
M87. It is interesting to note that the over-density required for a blob to be unstable is 
moderate: 5c ~ 0.8 — 3 if the blobs have sizes a m 10 — 100 pc. From Figure 8 we also note 
that the blobs have a minimum stable size. As it is seen from Equation (14), the heating 
term is proportional to the fall velocity, which is v t oc a 1//2 , where a is the blob radius. On 
the other hand, the radiative cooling term does not depend on the blob size. Therefore, 
we expect that if a blob is too small (i.e. below a critical threshold size) the gravitational 
term is insufficient to heat the blob: the radiative cooling prevails altogether, and the blob's 
over-density increases monotonically with time. 

Before ending this Section, we must notice that for the very over-dense blobs the as- 
sumption of isobaric cooling may break down, since the sound waves inside the blob become 
too slow to pressurize it against the ambient gas. The evolution is now isochoric, and the 
blob density grows less than the amount predicted by our isobaric model (e.g. Burkert & 
Lin 2000). If this transition occurs before the blob has reached the inner density plateau, 
its density could be not high enough to avoid the thermal stabilization. In a more typical 
situation the transition to the isochoric regime occurs when the blob has already reached 
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a sizeable over-density, so this effect is generally not very important, and the qualitative 
conclusions drawn under the isobaric assumption are not altered; the quantitative results 
may be different by small amounts for the small blobs. 

3.1. The Issue of Angular Momentum 

In our qualitative sketch we have so far omitted any reference to the angular momentum 
of the infalling material. A too large angular momentum might prevent the flow from 
approaching the central black hole: the flux would merely stagnate, cool down and condense 
in filaments or stars. The AGN fueling is cut off altogether, which makes the feedback 
impossible (see Cowie et al. 1980, for a thorough analysis of the thermal instability of a high 
angular momentum flow). The existence of a circumnuclear disc around M87 shows that the 
the flux possesses an amount of angular momentum, so the question is whether this angular 
momentum is high enough to jeopardize the feedback. We argue it is not. 

First of all, as we will discuss below, the blobs are expected to form and accrete only in a 
region of the same extension as the inner gas density plateau (5 — 30 kpc). The distribution 
of angular momentum at larger radii is immaterial to the present discussion, and in the rest 
of this Section we only refer to the ICM within this central region. 

The non-linear perturbation spectrum may stem directly from ICM disturbances driven 
by an early AGN activity, but also from galaxies mass-stripping (Soker et al. 1991). Since 
the galaxies do not have an ordered bulk motion, the blobs stripped from them are also 
unlikely to organize in an ordered flow with high net angular momentum. Besides, if a blob 
from a galaxy is injected with a high angular momentum, it is likely to loose most of it on 
account of its friction with the surrounding ICM. Therefore, even if a circular flow like a disc 
may form, it cannot be very large, as the example of M87 shows. 

A flow with specific angular momentum I circularizes at the radius -R C i rc ~ 1 2 /GMbh, 
being Mbh the central black hole mass. We estimate -R C i rc as follows. Let us consider a 
cold blob of radius a initially at the distance R from the central black hole. The balance 
between gravitational attraction and the friction force quickly brings the blob to the terminal 
velocity v t given by Equation (6). Then, its orbital angular momentum is about I ~ v t R 
and Equation (6) gives -R C irc ~ a, 5. The circularization radius is expected of the same size 
as the non-linear density blobs, since 5 is expected to be of the order of a few. Rough as it 
is, this estimate is in fair agreement with the actual size of the circumnuclear disc of M87 
Rd ~ 10 2 pc (Harms et al. 1994; Ford et al. 1994, see also below). We note that, statistically, 
many of the cold blobs will start with very low angular momentum I <C v t R. These will be 



-14- 



accreted directly to the black hole vicinity. 

As an example we consider the CF cluster M87. First we note that M87 has a large 
Keplerian disc at its center. Harms et al. (1994) estimate the black hole mass to be Mbh ~ 
2.4 x 10 9 M Q from their observation of a disc with a radius of ~ 20 pc. Their optical HST 
image shows the disc to be ~ 3.5 times larger, i.e., ~ 70 pc (Ford et al. 1994). 

The cooler gas at the center of the Virgo cluster is at a temperature of ~ 1 keV, which 
implies for the above black hole mass a Bondi accretion radius of ~ 80 pc (Churazov et al. 
2002) 1 . A Bondi accretion radius as large as the disc around the black hole of M87 further 
suggests that the simple Bondi accretion flow (Churazov et al. 2002; Nulsen 2004) does not 
hold; the accreted material has a larger angular momentum, and may come from much larger 
radii. 

To further elaborate on the proposed model, we plot in Figure 2 the free fall time in 
M87, as well as some relevant cooling times. In that Figure the crosses mark the cooling 
time of the ambient gas at several radii; the asterisks, the filled and the empty triangles are 
the cooling times of blobs with local over-densities with respect to the ambient of Sp/p a = 1, 
dp/ Pa = 3 and Sp/p a = 10, respectively. Pressure equilibrium is assumed between the blobs 
and the ambient medium: since the blobs temperature is always above T ^ 0.1 keV, the 
sound waves are fast enough to ensure this equilibrium. The solid line represents the free 
fall time in the gravitational potential of M87; the dashed line is the fall time of a blob of 
radius a = 300 pc and over-density Sp/p a = 1. In this case, the fall velocity is given by 
the balance between the gravitational acceleration and the drag force. The actual fall time 
lies in between these time-scales (see the next Section). The blobs for which t coo i m tf a n 
have moderate over-densities Sp/p a = 1 — 2, and therefore temperatures of 1/2 — 1/3 of the 
ambient gas. This value is not far from the values found by Molendi (2002) for the cool 
component in M87, and in good agreement with the cool component of NGC 5044 (Buote 
et al. 2003). 

An important final remark is in order. As Nulsen (1986) pointed out, in the absence of 
a cohesive force a blob would be torn apart by the ram pressure in a characteristic time 

U = -( £ -) V2 - 10 7 ( -?—) ( I" r ) " ( A V2 yr, (17) 
vt\PaJ \100 pcj \10 kms-V \p a J K ' 

where a is the blob's radius and v t its terminal velocity, given by Equation (6). This may 
be considerably shorter than the time taken by the blob to fall to the center, and we must 



1 Macchetto et al. (1997) estimate a somewhat larger mass for the central black hole in M87, namely 
^bh ~ (3-2 ± 0.9) x 10 9 M Q . In this case, the Bondi accretion radius is slightly larger, about 100 pc. 
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therefore assume that some kind of cohesive force (like a magnetic tension) is at work to 
prevent the blob disruption. In any case, the proposed scenario works also for the smaller 
blobs forming from the fragmentation of a larger blob. 

4. Summary 

This paper deals with heating the intra-cluster medium (ICM) in galactic and cluster 
cooling flows (CF) by an active galactic nucleus (AGN) sitting at the cooling flow center. 
As was shown by many papers, the heating is most likely to take place via a feedback 
mechanism, where the ICM cooling enhances the AGN activity, which in turn heats the 
ICM and quenches the cooling flow. 

Most previous papers (see § 1), assume that the central black hole accretes mass only 
from the ICM in its immediate neighborhood, basically via a Bondi-like accretion flow. In 
these models the feedback occurs as the ICM cools to a temperature of about 1 keV, and 
the ICM does not need to cool to low temperatures. We term these hot-feedback models. We 
examined three papers based on hot-feedback heating (Omma & Binney 2004; Ruszkowski 
& Begelman 2002; Hoeft & Briiggen 2004). We argued that the models worked out in these 
papers do not fit the general temperature profiles of CFs (Figure 1), and/or require fine 
tuning. We further argued that more generally, in the Bondi-type accretion flow of hot gas, 
the accretion rate is determined mainly by the conditions very close to the central black 
hole, and that this may result in unstable cooling of the regions further out. 

We therefore proposed (§ 2) that the feedback occurs with the entire cool inner region, 
r < 5 — 30 kpc, in what we term a cold-feedback model. In the proposed scenario non-linear 
over-dense (5p/p a > 2, or p/p a > 3) blobs of gas cool fast and are removed from the ICM 
before the next major AGN heating event in their region. It is important to note that an 
AGN burst can take place and heat other regions, since the jets and/or bubbles may expand 
in other directions as well. The typical interval between such heating events at a specific 
region is ~ 10 8 yr. Some of these blobs cool and sink toward the central black hole, while 
others may form stars and cold molecular clouds. 

Four conditions should be met in the inner region participating in the feedback heating. 

1. In order for the blob not to reach a point where its density equals the ambient density 
as it sinks, the ICM density profile should not be too steep. This implies that the 
relevant dense blobs form in the cluster core, where the density profile is shallow. In 
the quantitative example used here for M87 this is the region r < 5 kpc, while in 
more typical clusters it is larger, e.g., r < 30 kpc in A2052. We note that the lower 
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segment of magnetic flux loops can be prevented from reaching the stabilizing point 
by the upward force of the magnetic tension inside the loop (Soker 2004). Therefore, 
some perturbations can be formed at large distances, where density profile is steep, 
and still cool to low temperature and feed the central black hole. 

2. Non-linear perturbations are required. These presumably formed mainly by previous 
AGN activity, e.g. jets and radio lobes. 

3. The cooling rate of these non-linear perturbations is short relative to few times the 
typical interval between successive AGN outbursts. 

4. The blobs must not be evaporated by thermal conduction before they are delivered 
to the AGN. This requires a strong suppression of thermal conduction, which may be 
done by the magnetic fields observed in the radio lobes in several cooling flow clusters. 

The first and the third condition, which are not completely independent of each other, require 
that the initial ICM cools by a factor of a few before the feedback starts operating, and the 
second condition requires that the inner region must be disturbed. 

Finally, in § 3 we have calculated the falling time and cooling time of dense blobs. The 
results have then been applied to the cooling-flow cluster M87. 

The cold-feedback model has the following implications and predictions (for more details 
and references see § 1). 

1. In the cold-feedback models the optical filaments observed in many CF-clusters and 
the cooler molecular gas detected via CO observations come from cooling ICM (with 
some amount possibly from stripping from galaxies). 

2. In cold feedback models, some X-ray emission from gas at temperatures < 10 7 K is 
predicted to exist, much more than in many other AGN heating models, but at a level 
more than an order of magnitude below that in old versions of the CF model, but 
compatible with the moderate CF model. We stress that in the cold-feedback heating, 
cooling flows do exist. Such gas cooling to below X-ray emitting temperatures was 
found recently in two CF clusters (Abell 2597: Morris & Fabian 2005; Abell 2029: 
Clarke et al. 2004). 

3. The feeding of the central black hole with cold gas in the cold feedback models makes 
the process similar in some aspects to that of AGN in spiral galaxies. Therefore, the 
outflow can be similar (Soker & Pizzolato 2005). 
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4. It is possible that in the cold feedback model a substantial fraction of gas that cooled 
to low temperatures and was accreted to the accretion disc around the central black 
hole, is injected back to the ICM at non-relativistic velocities (Soker & Pizzolato 2005). 
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Fig. 1. — The temperature profiles for a sample of Clusters observed with Chandra (C) 
or XMM-Newton (X): A478 (Sun et al. 2003b, C), A496 (Tamura et al. 2001a, X), A1068 
(Wise et al. 2004, C), A1795 (Ettori et al. 2002, C), A1835 (Schmidt et al. 2001, C), A1991 
(Sharma et al. 2004, C), A2029 (Lewis et al. 2003, C), A2052 (Blanton et al. 2001, C), A2199 
(Johnstone et al. 2002, C), A4059 (Choi et al. 2004, C), M87 (Ghizzardi et al. 2004, C+X), 
Perseus (Schmidt et al. 2002, C), and finally the groups NGC 1550 (Sun et al. 2003a, C) and 
NGC 5044 (Buote et al. 2003, C). For the sake of readability, the clusters with an approxima- 
tively power-law temperature profiles and those with a central temperature floor have been 
plotted in two different panels (upper and lower, respectively). Where necessary, the radii 
from the original papers have been corrected for a cosmology with H = 70 km/s/Mpc. 
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Fig. 2. — A comparison between the free fall time, the terminal- velocity fall time and the 
isobaric cooling times for M87. The free fall time scale tg ~ (2r/g) 1//2 is shown by the 
solid line; the terminal- velocity fall times t t ~ r/v t refer to a blob of radius a = 300 pc 
and an over-density Sp/p a = 1 (dashed line) and Sp/p a = 3 (dotted line). The gravitational 
acceleration has been calculated from the deprojected density and temperature profiles, with 
the additional hypothesis of hydrostatic equilibrium. For the data and the procedure, see 
Ghizzardi et al. (2004), and references therein. The cooling times have been calculated with 
the deprojected data of Ghizzardi et al. (2004), assuming an average metal abundance of 
Z/Zq = 1.0 (expressed in Anders & Grevesse 1989, solar units), with the cooling function 
provided by Sutherland & Dopita (1993). The crosses refer to the ambient gas, asterisks to 
gas blobs with local over-densities Sp/p a = 1, the filled and empty triangles respectively to 
gas blobs with local over-densities Sp/p a = 3 and Sp/p a = 10 with respect to the ambient 
gas. Pressure equilibrium between the blobs and the ambient is assumed. The plot does not 
extend beyond 0.4 kpc because the density and temperature profiles we have taken from 
Ghizzardi et al. (2004) and that we used to calculate the characteristic time scales do not 
push to smaller radii. 
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Fig. 3. — The evolution of blobs with ao = 100 pc, tq = 20 kpc, but different initial over- 
densities 5 = (p — Pa)/ Pa' = 1-5 (solid line), 5 = 2.0 (dashed line), 5 = 2.5 (dot-dashed 
line), 5q = 3.0 (dotted line). Below a critical threshold 5c of the initial over-density 5q (in 
this case 2 < 5c < 2.5), the blobs are stabilized, above this limit they condense in a short 
(cooling) time. Here and in the following Figures the ICM properties reproduce those of 
M87, whose data have been derived from Ghizzardi et al. (2004). 
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Fig. 4. — The distance from the cluster center (upper panel) and the fall velocity (lower 
panel) of two blobs with initial over-densities 5 = 2.0 (solid line) and 5 = 2.5 (dashed line). 
Both the blobs have the initial radius ao = 100 pc, and have been released from rest at the 
distance tq = 20 kpc from the cluster center. 
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Fig. 5. — The temperature evolution of two blobs with the same characteristics as in Figure 4. 
Also the line styles are the same as in Figure 4. 
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Fig. 6. — The evolution of the over-density as a function of the release radius: r = 2 kpc 
(solid line), r = 5 kpc (dashed line), r = 10 kpc (dash-dotted line), and r = 20 kpc (dotted 
line). In all cases the blob has an initial over-density 5q = 2, and an initial radius Qq = 100 pc. 
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Fig. 7. — The total number particle density (upper panel) and the temperature profile (lower 
panel) of the ambient gas in M87. These profiles have been been borrowed from Ghizzardi 
et al. (2004). 
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Fig. 8. — The critical density 5c as a function of the blob size a for different initial positions 
of the blob: r = 2 kpc (solid line), r = 5 kpc (dashed line), r = 10 kpc (dash-dotted 
line), and r = 40 kpc (dotted line). The blobs on the left part of the plot are unstable, 
and evolve to larger values of 5p/ p a - The blobs on the right part are stable, and evolve 
to smaller over-densities. The blobs with sizes in the range a m 10 — 100 pc are unstable 
provided that their over-densities are in the moderately non-linear regime Sp/p a ~ 0.8 — 3, 
or the that they have been released from very small radii. For consistence with our neglect 
of the cooling of the ambient gas, we should consider only the part of the curves above the 
horizontal dashed line, corresponding to a blob cooling time of about 30% of the ambient 
gas. Blobs with over-densities below this (rather arbitrary) limit have cooling times closer 
to the ambient gas. The 30% confidence limit has been calculated by comparing the cooling 
time t a oc T a /n a A(T a ) for the ambient gas with the corresponding expression for the blob 
cooling time tb- We assume pressure equilibrium between the ambient and the blob, so 
T a = (1 + 5) Tb, where Tb is the blob temperature. By approximating the cooling function 
with a power law A oc T a , we obtain t a /t b = (1 + 5) 2 ~ a . Our approximation t b <C t a requires 
5 3> 1. If we demand tb < 30% t a and take a ~ in the temperature range considered here, 
we obtain the plotted confidence limit 5 > 0.8. 




